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ABSTRACT 

We have used observations of Cygnus X-l from the Compton Gamma-Ray Observatory 
(CGRO) and BeppoSAX to study the variation in the MeV 7 -ray emission between the hard 
and soft spectral states, using spectra that cover the energy range 20 keV up to 10 MeV. These 
data provide evidence for significant spectral variability at energies above 1 MeV. In particular, 
whereas the hard X-ray flux decreases during the soft state, the flux at energies above 1 MeV 
increases, resulting in a significantly harder 7 -ray spectrum at energies above 1 MeV. This be¬ 
havior is consistent with the general picture of galactic black hole candidates having two distinct 
spectral forms at soft 7 -ray energies. These data extend this picture, for the first time, to energies 
above 1 MeV. We have used two different hybrid thermal/non-thermal Comptonization models 
to fit broad band spectral data obtained in both the hard and soft spectral states. These fits 
provide a quantitative estimate of the electron distribution and allow us to probe the physical 
changes that take place during transitions between the low and high X-ray states. We find that 
there is a significant increase (by a factor of ~ 4) in the bolometric luminosity as the source 
moves from the hard state to the soft state. Furthermore, the presence of a non-thermal tail in 
the Comptonizing electron distribution provides significant constraints on the magnetic field in 
the source region. 

Subject headings: accretion, accretion disks — black hole physics — gamma rays: observations — stars: 
individual (Cygnus X-l) — X-rays: stars 


1. Introduction 

High energy emission from galactic black hole 
candidates (GBHCs) is characterized by variabil- 
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ity on time scales ranging from msec to months. 
In the case of Cygnus X-l, it has long been rec¬ 
ognized that, on time scales of several weeks, the 
soft X-ray emission (~ 10 keV) generally varies be¬ 
tween two discrete levels (e.g., Priedhorsky et al. 
1983; Ling et al. 1983; Liang & Nolan 1983). The 
source seems to spend most (~ 90%) of its time 
in the so-called low X-ray state , characterized by 
a relatively low flux of soft X-rays and a relatively 
high flux of hard X-rays (~ 100 keV). This state 
is sometimes referred to as the hard state, based 
on the nature of its soft X-ray spectrum. On occa¬ 
sion, it moves into the so-called high X-ray state, 
characterized by a relatively high soft X-ray flux 
and a relatively low hard X-ray flux. This state 
is sometimes referred to as the soft state, based 
on the nature of its soft X-ray spectrum. There 
are, however, some exceptions to this general be- 
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havior. For example, HEAO-3 observed, in 1979, 
a relatively low hard X-ray flux coexisting with a 
low level of soft X-ray flux (Ling et al. 1983, 1987). 
Ubertini et al. (1991) observed a similar behavior 
in 1987. 

Observations by the BATSE, OSSE, COMP- 
TEL and EGRET instruments on the Compton 
Gamma-Ray Observatory ( CGRO ), coupled with 
observations by other high-energy experiments 
(e.g., SIGMA, ASCA and RXTE) have provided 
a wealth of new information regarding the emis¬ 
sion properties of galactic black hole candidates. 
One important aspect of these high energy radia¬ 
tions is spectral variability, observations of which 
can provide constraints on models which seek to 
describe the global emission processes. Based on 
observations by OSSE of seven transient galactic 
black hole candidates at soft y-ray energies (i.e., 
below 1 MeV), two y-ray spectral shapes have 
been identified that appear to be well-correlated 
with the soft X-ray state (Grove et al. 1997, 1998; 
Grove 1999). In particular, these observations de¬ 
fine a breaking y-ray spectrum that corresponds 
to the hard (low) X-ray state and a power-law y- 
ray spectrum that corresponds to the soft (high) 
X-ray state. 

A thorough understanding of the nature of 
these systems requires modeling that cannot only 
explain the individual spectra, but can also ex¬ 
plain the transitions between the various spectral 
states (e.g., Grove et al. 1998; Liang 1998; Pouta- 
nen 1998a). In recent years, a general theoretical 
picture of the accretion flow in Cygnus X-l has 
emerged which appears to provide a reasonable 
explanation of the spectral data in both the low 
and high X-ray states. This model includes an in¬ 
ner optically-thin, geometrically-thick advection- 
dominated accretion flow (ADAF) surrounded by 
an outer, geometrically-thin, optically-thick accre¬ 
tion disk (Esin et al. 1998). The outer disk is 
characterized by a blackbody spectrum. ADAF 
flows (e.g., Narayan 1996) are characterized by 
their relatively low radiative efficiencies and by a 
two-temperature structure, with the ions nearly 
virial at T) ~ 10 12 K and the electrons at T e ~ 10 9 
K. The high temperature of the ADAF leads to an 
extended, quasi-spherical geometry. Hot optically- 
thin ADAFs exist only below a certain critical ac¬ 
cretion rate. The transition radius between the 
ADAF and the thin disk therefore depends on the 


accretion rate. At higher accretion rates, where it 
is more difficult to support the ADAF, the tran¬ 
sition radius moves to smaller radii, closer to the 
black hole. The ADAF region is largely responsi¬ 
ble for the hard X-ray flux (~ 20 — 100 keV), while 
the outer thin disk is generally responsible for the 
soft X-ray flux (~ 2 — 10 keV). 

In the context of this general model, the spec¬ 
tral state of Cygnus X-l depends on the ac¬ 
cretion rate. At low accretion rates, the inner 
ADAF extends out to a transition radius of ~ 100 
Schwarzschild radii (Esin et al. 1998). In this 
configuration, the ADAF region makes a signif¬ 
icant contribution to the hard X-ray flux. At 
higher accretion rates, it becomes more difficult 
to support the ADAF. The ADAF region therefore 
shrinks, and the transition radius moves inward, 
although there may exist a low-density ADAF 
corona surrounding the thin disk (Narayan et al. 

1998) . The level of hard X-ray flux decreases due 
to the smaller volume of the ADAF region, while 
the level of soft X-ray flux increases due to the 
larger size of the thin disk region. In this sce¬ 
nario, the hard state corresponds to a relative low 
accretion rate, with the spectrum dominated by 
the ADAF region, and the soft state corresponds 
to a relatively high accretion rate, with the spec¬ 
trum dominated by the blackbody of the outer 
thin disk region. Small changes in the accretion 
rate (on the order of 10-15%) may be sufficient 
to trigger a transition between the hard and soft 
states (Esin et al. 1998). 

The ADAF model described above provides a 
consistent framework for understanding the essen¬ 
tial dynamics and spectra of black hole accretion 
flows. In the context of this framework, how¬ 
ever, simple thermal Comptonization models ap¬ 
pear unable to account for all of the spectral fea¬ 
tures, especially the hard power-law tail that is 
seen at energies above ~ 600 keV (Gierlinski et al. 

1999) . Poutanen & Coppi (1998) used a geometry 
similar to that described above (Poutanen et al. 
1997) and assumed some (unspecified) source of 
non-thermal electrons that remains constant dur¬ 
ing the spectral state transitions. This suggests 
that the non-thermal component may play a more 
significant role, especially at higher energies, dur¬ 
ing the high X-ray state, where the ADAF contri¬ 
bution is suppressed. 

Hybrid thermal/non-thermal plasmas have of- 


2 



ten been successfully used to model the observed 
data (e.g., Gierlinski et al. 1999; Poutanen & 
Coppi 1998). Based on the assumption that 
the spectrum results from inverse Compton scat¬ 
tering of a thermal photon spectrum by ener¬ 
getic electrons, the underlying electron population 
could be described as a combination of a ther¬ 
mal Maxwellian and a power-law tail extending 
to higher energies. The presence of a non-thermal 
component is often assumed a priori, without any 
specific model to explain the origin, although the 
existence of such distributions is clearly estab¬ 
lished in the case of solar flares (e.g., Coppi 1999) 
and it is therefore natural to expect that similar 
distributions exist elsewhere in the universe (e.g., 
Crider et al. 1997; Gierlinski et al. 1997; Poutanen 
& Svensson 1996; Poutanen 1998a; Poutanen & 
Coppi 1998; Coppi 1999). Others have considered 
physical mechanisms by which non-thermal elec¬ 
tron distributions might be developed. For exam¬ 
ple, both stochastic particle acceleration (Derrner 
et al. 1996; Li et al. 1996) and MHD turbulence 
(Li & Miller 1997) have been proposed as mecha¬ 
nisms for directly accelerating the electrons. The 
ion population might also contribute to the non- 
thermal electron distribution in the case where a 
two-temperature plasma develops (e.g., Dahlbacka 
et al. 1974; Shapiro et al. 1976; Chakrabarti & 
Titarchuk 1995). With ion population tempera¬ 
tures approaching kTi ~ 10 12 K, n° production 
from proton-proton interactions may take place 
(e.g., Eilik 1980; Eilik & Kafatos 1983; Mahade- 
van et al. 1997). The 7r° component may then 
lead, via photon-photon interactions between the 
7r°-decay photons and the X-ray photons, to pro¬ 
duction of energetic (nonthermal) e + — pairs. 
Jordain & Roques (1994) used this concept to fit 
the hard X-ray tails of not only Cygnus X-l, but 
also GRO J0422+32 and GX 339-4, as measured 
by both SIGMA and OSSE. While retaining a 
standard thermal Comptonization spectrum (Sun- 
yaev & Titarchuk 1980) to explain the emission at 
energies below 200 keV, they used tt° production 
to generate the nonthermal pairs needed to fit the 
spectrum at energies above ~ 200 keV. 

The power-law spectra seen in the high X-ray 
state have also been modeled as resulting from 
bulk-motion Comptonization (e.g., Ebisawa et al. 
1996b; Titarchuk et al. 1997; Laurent & Titarchuk 
1999). In this model, the flow becomes quasi- 


spherical within the inner-most stable orbit. The 
nearly relativistic flow of the free-falling electrons 
gives rise to the Comptonization of ambient pho¬ 
tons. This model predicts power-law spectra, with 
a slope that depends on the mass acretion rate. 
The difficulty with this model is that it predicts 
spectral sharp cutoffs below 500 keV, a result that 
is clearly inconsistent with the observed spectra. 
Although we cannot rule out bulk motion Comp¬ 
tonization as a contributor to the spectrum at 
lower energies, it is clearly not capable of account¬ 
ing for the high energy emission. 

Improvements in the theoretical modeling of 
spectral state transitions can be expected to arise 
from improved observations at energies above 600 
keV. It will be important to understand how this 
part of the spectrum, most likely dominated by 
non-thermal emission, changes during the spectral 
transition. Of particular interest will be observa¬ 
tions that can discern a clear cutoff in the spectra 
at high energies. The precise energy of the cutoff is 
a function of the compactness of the source region, 
since it is influenced by 7—7 opacity. A measure of 
the cutoff energy, possibly coupled with measure¬ 
ments of the 511 keV e ± annihilation line, will help 
constrain the compactness of the region responsi¬ 
ble for the emission and determine the extent to 
which e ± pairs may play a role in the emission 
region (Poutanen 1998a). 

Using hard state data collected during the first 
three years of the CGRO mission, McConnell et 
al. (2000a) compiled a broad-band hard state spec¬ 
trum of Cygnus X-l using contemporaneous data 
from all four instruments on CGRO (BATSE, 
OSSE, COMPTEL and EGRET). Unlike previ¬ 
ous broad-band studies, these data provided a 
measurement of the spectrum at energies above 
1 MeV. The resulting spectrum showed evidence 
for significant levels of non-thermal emission at 
energies out to 5 MeV. The spectral shape, al¬ 
though consistent with the so-called breaking 
spectral state (Grove et al. 1997, 1998) of the 
7-ray emission, was clearly not consistent with 
standard Comptonization models. The hybrid 
thermal/non-thermal model of Poutanen & Svens¬ 
son (1996) was used to fit the hard state data, with 
fits that indicated a thermal electron population 
with a temperature of ~ 90 keV and a high energy 
power-law electron component with a spectral in¬ 
dex of ~ 4.5. 
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In May of 1996, a transition of Cyg X-l into a 
soft state was observed by RXTE, beginning on 
May 10 (Cui 1997). The 2-12 keV flux reached a 
level of 2 Crab on May 19, four times higher than 
its normal value. Meanwhile, at hard X-ray ener¬ 
gies (20-200 keV), BATSE measured a significant 
decrease in flux (Zhang et al. 1997). Motivated by 
these dramatic changes, a target-of-opportunity 
(ToO) for CGRO, with observations by OSSE and 
COMPTEL began on June 14 (CGRO viewing pe¬ 
riod 522.5). Here we report on the results from an 
analysis of the CGRO data from this ToO obser¬ 
vation, incorporating the high energy results from 
COMPTEL. This includes a comparison with re¬ 
sults obtained from an updated analysis of CGRO 
soft state data, making use of the same data stud¬ 
ied previously by McConnell et al. (2000a). In § 2 
we describe the CGRO observations of Cygnus X- 
1 in its hard state. The data analysis is described 
in § 3, followed by a discussion of those results in 
§ 4 . 

2. Observations 

During its nine-year lifetime (1991-2000), the 
instruments on CGRO obtained numerous obser¬ 
vations of the Cygnus region. The COMPTEL 
experiment (Schonfelder et al. 1993), imaging the 
energy range from about 750 keV up to 30 MeV, 
collected an extensive set of data, in part due to 
its rather large (~ ir steradian) field-of-view. The 
COMPTEL data currently provide the best avail¬ 
able source of data for studies of Cygnus X-l at 
energies above 1 MeV. 

The 20-100 keV time history of Cygnus X-l, as 
derived from BATSE occultation data, is shown 
in the center panel of Figure 1. The top panel 
of Figure 1 shows the 20-100 keV power-law spec¬ 
tral index, as derived from the BATSE occultation 
data. These data cover most of the CGRO mis¬ 
sion, from the launch in April of 1991 until the 
end of 1999. During the first few months of the 
CGRO mission (up until October of 1991), all¬ 
sky monitoring data from Ginga (1-20 keV) was 
available, showing that the source was in its low 
X-ray state during this period (Kitamoto et al. 
2000). From October of 1991 until December of 
1995, there were only sporadic pointed X-ray ob¬ 
servations of the soft X-ray flux from Cygnus X-l. 
It was not until the launch of RXTE, in Decem¬ 


ber of 1995, that continuous data on the soft X-ray 
flux once again became available. The data from 
the RXTE All-Sky Monitor (ASM) are shown in 
the lower panel of Figure 1, in the form of the 2-10 
keV count rate. 

The data shown in Figure 1 dramatically 
demonstrate the general X-ray behavior of Cygnus 
X-l. During the CGRO mission, Cygnus X-l 
spent about 90% of its time in the hard state. 
In this state, the soft X-ray flux (2-10 keV) is 
relatively low, while the hard X-ray flux (20-100 
keV) is relatively high. The spectral shape in 
the 20-100 keV energy band is a relatively hard 
power-law spectrum with a photon spectral index, 
r, near 1.8. The soft state was clearly observed 
during the CGRO mission on only two occassions. 
In each case, the soft state period lasted about 5 
months. The soft state is characterized by a rela¬ 
tively high level of soft X-rays (2-10 keV), a rela¬ 
tively low level of hard X-rays (20-100 keV) and a 
relatively soft spectrum in the 20-100 keV energy 
band (photon spectral index ~ 2.5). The soft state 
was first observed by CGRO in January of 1994, 
at a time (prior to the launch of RXTE) when 
there was no soft X-ray monitoring data available. 
(This transition is clearly seen in Figure 1 near 
TJD 9400.) A CGRO target-of-opportunity was 
declared (CGRO viewing period 318.1) so that all 
four CGRO instruments (not just BATSE) could 
collect data. Observations by COMPTEL showed 
no detectable level of emission. This null result, 
however, was consistent with an extrapolation of 
the E~ 2,7 power-law spectrum measured at hard 
X-ray energies by both BATSE (Ling et al. 1997) 
and OSSE (Phlips et al. 1996). 

The second observation of a soft state took 
place in May of 1996. The transition was first 
observed by RXTE, beginning on May 10 (Cui 
1997). The 2-12 keV flux reached a level of 2 
Crab on May 19, four times higher than its nor¬ 
mal value. Meanwhile, at hard X-ray energies (20- 
200 keV), BATSE measured a significant decrease 
in flux (Zhang et al. 1997). Motivated by these 
dramatic changes, a ToO for CGRO was declared 
and observations by OSSE and COMPTEL began 
on June 14 (CGRO viewing period 522.5). (Un¬ 
fortunately, the EGRET experiment was turned 
off during this viewing period, as part of an ef¬ 
fort to conserve its supply of spark chamber gas.) 
During the ToO, COMPTEL collected 11 days of 
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data (from June 14 to June 25) at a favorable as¬ 
pect angle of 5.3°. The X-ray flux time histories 
near the time interval associated with VP 522.5 
are shown in Figure 2. 

An early preliminary analysis of COMPTEL 
data from this second high-state observation re¬ 
vealed some unusual characteristics (McConnell et 
al. 2000b). The 1-3 MeV image (Figure 3) showed 
an unusually strong signal from Cygnus X-l when 
compared with other observations of similar ex¬ 
posure. The flux level was significantly higher 
than the average flux seen from earlier observa¬ 
tions (McConnell et al. 1994, 2000a). In the 1-3 
MeV energy band, the flux had increased by a fac¬ 
tor of 2.5, from 8.6(±2.7) x 10 -5 cm -2 s -1 MeV -1 
to 2.2(±0.4) x 10" 4 cm" 2 s" 1 MeV" 1 . The ob¬ 
served change in flux is significant at a level of 
2.6(7. In addition, unlike in previous measure¬ 
ments, there was no evidence for any emission at 
energies below 1 MeV. This fact is explained, in 
part, by a slowly degrading sensitivity of COMP¬ 
TEL at energies below 1 MeV due to increasing en¬ 
ergy thresholds in the lower (D2) detection plane. 
Part of the explanation, however, appears to be 
the presence of a much harder source spectrum. 
OSSE data collected during this period showed a 
photon spectrum similar to that observed in 1994 
(i.e., a power-law with an index of ~ 2.5), but at a 
higher intensity level, about a factor of two higher 
in overall normalization. The extrapolation of this 
more intense power-law spectrum is entirely con¬ 
sistent with the positive detection by COMPTEL. 

3. Data Analysis 

A more detailed description of the COMP¬ 
TEL data analysis is given in McConnell et al. 
(2000a). Here we provide only a brief overview. 
The COMPTEL image shown in Figure 3 is a 
maximum likelihood map derived from VP 522.5 
data integrated over the energy loss range of 1-3 
MeV. The contours represent constant values of 
the quantity —2 In A, where A is the likelihood ra¬ 
tio. In a search for single point sources, —2 In A has 
a chi-square distribution with 3 degrees of free¬ 
dom. (For instance, a 3(7 detection corresponds 
to —2 In A = 14.2.) Cygnus X-l is clearly visible. 
The likelihood reaches a value of —2 In A = 30.1 
at the position of Cygnus X-l, which corresponds 
to a detection significance of 5.5o\ These same 


data were used to derive the lex, 2cr and 3a loca¬ 
tion confidence contours shown in Figure 4, which 
demonstrate the ability of COMPTEL to locate 
the source of emission. In defining constraints on 
the source location, —2 In A has a chi-square dis¬ 
tribution with 2 degrees of freedom. So the lex, 2cr 
and 3cr location confidence contours correspond to 
a change in likelihood of 2.3, 6.2, and 11.8, re¬ 
spectively. The contours reflect only the statis¬ 
tical uncertainties; systematic effects are not in¬ 
cluded. The COMPTEL flux results for VP 522.5 
are shown in Table 1. 

The analysis of COMPTEL data for a weak 
source (such as Cygnus X-l) involves generating 
an image for each of several energy bands, de¬ 
riving the source flux in each energy band, and 
subsequently combining these results into a spec¬ 
trum (e.g., McConnell et al. 2000a). The image 
generation process, in turn, requires an instru¬ 
ment point-spread-function (PSF) that is depen¬ 
dent on some assumed form for the incident pho¬ 
ton spectrum. Because the spectrum extraction 
relies on an assumed source spectrum (we typi¬ 
cally assume an E~ 2 power-law spectrum), it is 
not possible to analyze the COMPTEL data us¬ 
ing a simple response function to relate measured 
energy-loss count rates to the incident photon flux. 
We have therefore resorted to spectral fitting of 
COMPTEL data in photon space. We have pre¬ 
viously shown (McConnell et al. 2000a) that this 
approach to COMPTEL spectral analysis works 
fine for the range of spectra considered for Cygnus 
X-l. In other words, for the range of parameters 
considered here, there is no evidence of any signif¬ 
icant level of spectral compliance in the COMP¬ 
TEL spectral analysis. 

In previous work (McConnell et al. 2000a), we 
analyzed a contemporaneous set of CGRO data 
corresponding to the hard state of Cygnus X- 
1. That analysis was performed entirely in pho¬ 
ton space, using one deconvolved spectrum each 
for BATSE, OSSE and COMPTEL. The BATSE 
spectrum had been generated using the JPL En¬ 
hanced BATSE Occultation Package (EBOP; Ling 
et al. 1996, 2000), while the OSSE spectrum was 
based on the results of Phlips et al. (1996). The 
analysis was performed within XSPEC (Arnaud 
1996) to take advantage of the XSPEC analysis 
tools. In this case, however, the spectral data (pha 
files) were generated in units of photons cm' 2 s _1 
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and the response function matrices (rsp files) were 
generated as unit matrices. With these data, the 
spectral fits were effectively being performed in 
photon space. This approach greatly simplified 
the analysis effort. 

The analysis employed here for the soft state 
data from VP 522.5 represents a significant im¬ 
provement over that performed previously in gen¬ 
erating a broad-band gamma-ray spectrum for 
the hard state. Although the limitations of the 
COMPTEL data analysis remain, we have uti¬ 
lized more complete spectral response information 
for both BATSE and OSSE, using proper XSPEC 
pha and rsp files. The fundamental nature of the 
COMPTEL data (in particular, its reliance on an 
assumed PSF for extracting source counts) still 
precludes a proper XSPEC analysis of the COMP¬ 
TEL data. 

Since the BATSE EBOP processing has not 
been carried out for data collected after 1994, 
EBOP data for VP 522.5 is not available. Instead, 
we have used data derived from the BATSE team’s 
standard Earth occultation analysis (Harmon et 
al. 2002). The final BATSE spectrum represents 
a weighted average of the four forward-facing de¬ 
tectors. 

For consistency, in order to make a more useful 
comparison with the soft state data, we have re¬ 
peated the earlier hard state analysis (McConnell 
et al. 2000a) following the same procedures as we 
have used here for the VP 522.5 data. In particu¬ 
lar, the updated hard state analysis now also used 
BATSE data derived from the standard Earth oc¬ 
cultation technique (Harmon et al. 2002). Data 
from nine separate CGRO viewing periods were 
used (see Table 1 in McConnell et al. 2000a). In 
each case a weighted average spectrum was derived 
from the data for use in the final analysis. 

Finally, our most recent analysis incorporates 
the BATSE data down to 20 keV and OSSE 
data down to 50 keV. Previously, we had used 
only those data above 200 keV. The lower en¬ 
ergy threshold of this analysis improves the abil¬ 
ity of our fits to constrain the spectral models. 
At the same time, the lower threshold may also 
make the analysis more sensitive to systematic 
uncertainties in the low energy response of both 
BATSE and OSSE. The OSSE (Johnson et al. 
1993) data include energy-dependent systematic 
errors (estimated from the uncertainties in the 


low-energy calibration and response of the detec¬ 
tors using both in-orbit and prelaunch calibration 
data), which are most important at the lowest en¬ 
ergies, ~ 3% at 50 keV, decreasing to ~ 0.3% at 
> 150 keV. To the BATSE data, we added a 5% 
systematic error. The COMPTEL data have rela¬ 
tively large statistical errors, and thus no system¬ 
atic error was added. 

3.1. The Average Hard (Low) State Spec¬ 
trum 

X-rays from Cygnus X-l in the hard state 
are well modeled by thermal Comptonization and 
Compton reflection (Gierlinski et al. 1997; Di 
Salvo et al. 2001; Frontera et al. 2001). Thus, 
we fit the joint data by this model, but also al¬ 
lowing for a tail at high electron energies rather 
than the Maxwellian cutoff. This is similar to the 
approach used in McConnell et al. (2000a), except 
that here we have added the Compton reflection 
component, which is important at energies below 
a few hundred keV. 

We first fit the data using the Comptonization 
model (compps) 1 of Poutanen & Svensson (1996), 
assuming a spherical source geometry. The elec¬ 
trons have the total Thomson optical depth of 
r. Their distribution in this model is Maxwellian 
with an electron temperature, kT, up to a Lorentz 
factor, 7 m i n , above which it is a power law with an 
index, p. The power law extends to a large Lorentz 
factor, 7 max - The precise value of 7 max , however, 
has little effect on the fit to our data as long as 
7^ ax times the seed photon energy is > 10 MeV. 
Given that the seed photons peak at a fraction of 
keV (Ebisawa et al. 1996a; Di Salvo et al. 2001), 
we assume 7max = 10 3 . The Comptonization spec¬ 
trum is then Compton-reflected from a cold slab 
(presumably an accretion disk) subtending a solid 
angle, H (Magdziarz & Zdziarski 1995). A disk in¬ 
clination of i = 45° is assumed (as in Gierlinski et 
al. 1997; Frontera et al. 2001). For this model, we 
found H is not constrained by our data and we kept 
it fixed at a typical value Q/2ir = 0.5 (Gierlinski et 
al. 1997; Gilfanov et al. 1999; Di Salvo et al. 2001). 
This value also follows from our fit below using 
another theoretical model. The seed photons for 
Comptonization are assumed to be a (multicolor) 
blackbody emission of the disk with the maximum 

Available at ftp://ftp.astro.su.se/pub/juri/XSPEC/COMPPS 
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blackbody temperature of kT s = 0.2 keV, which 
approximately corresponds to a single blackbody 
temperature of ~ 0.13-0.15 keV obtained in the 
fits of Ebisawa et al. (1996a) and Di Salvo et al. 
( 2001 ). 

During the data analysis, we find some resid¬ 
ual discrepancies between the different data sets, 
as expected for different instruments. First, the 
BATSE spectrum has a slightly higher normaliza¬ 
tion than the OSSE one, and thus we allow their 
relative normalization to be free in the fits. Fur¬ 
thermore, the BATSE spectrum is systematically 
slightly softer, by AT ~ 0.1 (where V is the pho¬ 
ton power law index), than the OSSE spectrum. 
We find that multiplying the BATSE model by 
an additional power law with that AT leads to a 
reduction of y 2 /*' from 71/48 to 37/47, which is 
highly significant at the 2 x 10~ 8 level, using the 
F-test (Bevington & Robinson 1992). Thus, we 
apply this correction, fixing AT hereafter at the 
best-fit value. The best-fit ratio of the BATSE 
and OSSE fluxes at 100 keV is then 1.26 (but 
higher and lower at higher energies and lower en¬ 
ergies, respectively). Furthermore, we find that 
the COMPTEL data at ~ 1 MeV appear to have 
somewhat higher normalization than the BATSE 
and OSSE ones. The best-fit relative normaliza¬ 
tion is ~ 1.5, consistent with results for the soft 
state (see below). Given the limited statistics of 
the COMPTEL hard-state data, we fix that rela¬ 
tive normalization at 1.5. We note that this yields 
a conservative estimate of the amplitude of the 
nonthermal tail in Cygnus X-l. 

The compps fit results are given in Table 2. 
They are similar to the preliminary results of Mc¬ 
Connell et al. (2000a), who neglected Compton re¬ 
flection. We note that the fitted electron distribu¬ 
tion is allowed to be significantly different from a 
pure Maxwellian, with a power law tail beginning 
at a rather low energy. This reflects the fact that 
an arbitrary electron distribution peaked at low 
energies yields Comptonization spectra relatively 
similar to those from a pure Maxwellian (Ghis- 
cllini et al. 1993), apart from the high-energy tail 
in the former case. 

We then fit the same data using a different hy¬ 
brid Comptonization model, eqpair (Coppi 1992, 
1999; Poutanen & Coppi 1998; Gierlinski et al. 
1999). Unlike compps, which assumed the form 
of the steady-state electron distribution, eqpair 


calculates that distribution self-consistently as¬ 
suming instead the electron acceleration to be a 
power law with an index, r; n j between 7 m i n and 
7max- The accelaration takes place in a back¬ 
ground thermal plasma with a Thomson optical 
depth of ionization electrons, t\. The steady-state 
electron distribution consists then of a Maxwellian 
at the temperature, kT, calculated from the bal¬ 
ance of Compton and Coulomb gains and losses, 
and the optical depth, r, with r — r, due to the 
e± pair production. The nonthermal steady-state 
electron distribution is calculated from Coulomb 
and Compton losses of both the accelerated elec¬ 
trons and e ± pairs produced at nonthermal ener¬ 
gies. That distribution, in general, does not have 
a power law form. Unlike compps, the electron dis¬ 
tribution is now a sum of the Maxwellian and the 
nonthermal part (rather than being a Maxwellian 
up to 7min and then nonthermal). We assume 
7min = 1.5 and 7 max = 10 3 . Unlike the case of 
compps, where the value of 7 m i n determined the 
transition from the Maxwellian to the power law, 
that value has relative little effect on the fit now. 

The rates of microscopic processes per unit light 
travel time across the source depend in general 
on the plasma compactness, £ = Ca^ / (JZm e c 3 ), 
where £ is a power supplied to the hot plasma, 1Z 
is its characteristic size, and (Jt is the Thomson 
cross section (e.g., Svensson 1987). We then de¬ 
fine a hard compactness, £h, corresponding to the 
power supplied to the electrons, and a soft com¬ 
pactness, £ s , corresponding to the power in soft 
seed photons irradiating the plasma (which are as¬ 
sumed to be emitted by a blackbody disk). The 
compactnesses corresponding to the electron ac¬ 
celeration and to a direct plasma heating (i.e., in 
addition to Coulomb energy exchange with non¬ 
thermal and Compton heating) of the thermal 
e ± are denoted as £ n th and Ah, respectively, and 
fh = f'nth + £th- Details of the model are given in 
Gierlinski et al. (1999). 

The eqpair fit results are given in Table 2. Fig¬ 
ure 5 shows the spectrum. Figure 6 shows the 
model spectral components over the broad energy 
range from 0.1 keV to 100 MeV. The nonthermal 
high-energy tail starts at ~ 1 MeV. 

The bolometric flux derived from the eqpair 
model is very similar to that in the compps model. 
Some differences between the values of kT and r 
may be attributed to the different treatment of 
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the microphysics (see above). Also, the eqpair 
model gives a somewhat better fit to the data than 
the compps model. In addition, the eqpair model 
provides a better constraint on the value of Q,/2tt 
(Table 2). The better fit most likely reflects the 
fact that more physical processes are accounted 
for by eqpair than by compps. 

In particular, pair production is important at 
the best fit of the eqpair model, which accounts 
for r — r; > 0 in Table 2, and the associated injec¬ 
tion of nonthermal e* at low energies. The latter 
effect leads to a softening of the nonthermal spec¬ 
tra (Svensson 1987) and explains the relatively low 
value of Pi n j (without pair production and in the 
Thomson regime, p would be ~ Ti„j + 1, e.g., Blu- 
mentlial & Gould 1970). On the other hand, at the 
lowest 4 allowed by the data, ~ 0.2, the plasma 
compactness is so small that we find basically no 
pair production, i.e., t = t\. Thus, the present 
data do not resolve the issue of the role of pair 
production conclusively. 

The bolometric luminosity of the average hard- 
state spectrum, L, equals ~ 1% of the Eddington 
luminosity, Le — 1.5(M/M 0 ) x 10 38 erg s _1 , as¬ 
suming isotropy, a black-hole mass of M = 1OM 0 
and a distance of 2 kpc (see discussion and refer¬ 
ences in Gierlinski et al. 1999). The best-fit total 
compactness, -4+4 ~ 30, corresponds to the char¬ 
acteristic dimension of the plasma of 10 2 GM/c 2 
under the assumptions as above. The X-ray spec¬ 
trum is rather hard, with the amplification of the 
seed photons by the factor 4/4 — 17. Only a 
small fraction of the power supplied to the plasma, 
4th/ 4 — 0.08, is used for nonthermal electron ac¬ 
celeration. 

3.2. Broad-Band Spectrum in the Soft 
(High) State 

As in the case of the hard state, the broad¬ 
band spectrum of Cygnus X-l in the soft state is 
well fit by emission from a blackbody disk, Comp¬ 
ton scattering by thermal and nonthermal elec¬ 
tron components, and Compton reflection with the 
accompanying Fe Ka fluorescence line (Gierlinski 
et al. 1999; Frontera et al. 2001). However, un¬ 
like the hard state, the CGRO data alone (> 20 
keV) cannot determine the parameters of the ther¬ 
mal electron distribution. The reason for this is 
that, whereas scattering by the thermal electrons 
dominates up to several hundred keV in the hard 


state (see above, also Gierlinski et al. 1997), it 
dominates only up to ~ 10 keV in the soft state 
(Gierlinski et al. 1999, see below). Then, in the 
CGRO energy range, the spectrum is entirely due 
to the emission of the nonthermal electrons and 
Compton reflection. 

Thus, in order to determine the parameters 
of the electron distribution (including its thermal 
part) implied by the CGRO data in the soft state, 
we combine them with the BeppoSAX data from 
1996 June 22 (Frontera et al. 2001). These data 
cover roughly a 90 minute time span during the 
much longer 11-day CGRO observation (June 14- 
25). For the BeppoSAX observation, data from 
three instruments, LECS, HPGSPC and PDS, are 
usable (Frontera et al. 2001), extending the mea¬ 
sured energy range down to 0.5 keV. We allow for 
a free relative normalization of each set of spec¬ 
tral data with respect to that of OSSE. All the 
normalization factors are found to be ~ 1. 

We use the same two models (compps and 
eqpair) as for the hard state. However, since 
our data extend now down to ~ 0.5 keV, we let 
kT s free. For the same reason, we need to include 
the fluorescent Fe Ka emission, present in both 
states of Cygnus X-l. Since the line is produced 
by Compton reflection, we need to relate its flux 
to the strength of Compton reflection. We fol¬ 
low here results of George & Fabian (1991) and 
Zycki & Czerny (1994) and tie the line flux to 
ST/27r in such a way that the equivalent width with 
respect to the total continuum is ~ 120 eV when 
0/27r = 1. Both the line and the reflection contin¬ 
uum are assumed to come from an accretion disk 
extending down to 6 GM/c 2 (e.g., Gierlinski et al. 
1999) and with the reflection/fluorescence emissiv- 
ity following that of a standard thin disk (Shakura 
& Sunyaev 1973). This results in a relativistic 
smearing (Fabian et al. 1989) of both of those spec¬ 
tral components. The reflecting surface is allowed 
to be ionized (Gierlinski et al. 1999; Di Salvo et 
al. 2001), with the degree of ionization character¬ 
ized by the ionization parameter, £ = 4:TrFi on /n 
(where Fj on is the ionizing flux and n is the re¬ 
flector density), and at the temperature of ~ kT s . 
The elemental abundances are of Anders & Ebi- 
hara (1982). 

Both compps and eqpair models provide very 
good descriptions of our broad-band spectrum. 
Table 2 gives the fit results and Figure 7 shows the 



spectrum for the eqpair model. Figure 8 shows 
the spectral components of the eqpair fit to the 
spectrum. Both models predict the power-law-like 
emission extending with no cutoff up to 10 MeV, 
in agreement with the data. 

Strong Compton reflection with fl/2n ~ 1.3 is 
seen, similar to the results of Frontera et al. (2001) 
and those from RXTE (Gilfanov et al. 1999). A 
likely cause of Q > 2ir is relativistic anisotropy of 
Compton scattering (see a discussion in Gierlinski 
et al. 1999). Scattering by nonthermal electrons, 
forming a power-law like component, dominates 
a peaked component from thermal scattering at 
energies above several keV (as found by Gierlinski 
et al. 1999). 

Pair production is unimportant at the best fit. 
Thus, both models give results fully consistent 
with each other. The values of kT are virtually 
identical, and the small difference in the values of 
r is an artifact of the different treatment of the 
radiative transfer, and p ~ Fi„j + 1 (as expected 
for dominant Compton cooling in the Thomson 
regime, see §3.1). 

The bolometric luminosity is about 4 times that 
in the hard state, and it is ~ 0.04Le under the 
same assumptions as in §3.1. In contrast to the 
hard state, only a small fraction of the total lu¬ 
minosity, fh/(ih + 4) — 0.15, is emitted by the 
plasma outside the optically-thick accretion disk, 
although part of the disk emission is due to re¬ 
processing of the hard, plasma, emission (see dis¬ 
cussion of the energy balance in Gierlinski et al. 
1999). Also in contrast to the hard state, most, 
~ 0.7, of the power supplied to the plasma is used 
for nonthermal acceleration. Although the elec¬ 
tron temperature is very similar in both states, 
the optical depth in the soft state is <C that in the 
hard state. 

4. Discussion 

The COMPTEL data alone can be used to draw 
some important conclusions regarding the MeV 
variability of Cyg X-l. Most importantly, the flux 
measured by COMPTEL at energies above 1 MeV 
was observed to be higher (by a factor of 2.5) dur¬ 
ing the soft state (in May of 1996) than it was dur¬ 
ing the hard state (as averged over several CGRO 
observations). This is in contrast to the lower flux 
level observed at hard X-ray energies (i.e., near 


100 keV) during the soft state. The lack of any 
detectable emission by COMPTEL below 1 MeV 
(i.e., in the 750 keV to 1 MeV energy band) fur¬ 
ther suggests a hardening of the 7 -ray spectrum 
during the soft state. 

Inclusion of the BATSE and OSSE spectra adds 
considerably more information regarding the spec¬ 
tral variability. Whereas the low-state CGRO 
spectrum shows the breaking type spectrum that 
is typical of most high-energy observations of Cyg 
X-l (e.g., McConnell et al. 2000a), the high-state 
CGRO spectrum shows the power-law type spec¬ 
trum that is characteristic of black hole candidates 
in their high X-ray state. Our analysis of the soft 
state data from BATSE, OSSE and COMPTEL 
shows that the spectrum at these energies can be 
described by a single power-law with a best-fit 
photon spectral index of T = 2.58 ± 0.03. A sim¬ 
ilar spectrum had already been reported for this 
same time period (VP 522.5) based on indepen¬ 
dent studies with data from both BATSE (Zhang 
et al. 1997) and OSSE (Gierlinski et al. 1997, 
1999). A detailed study of the broadband soft 
state spectrum, based on data from ASCA, RXTE 
and CGRO/OSSE, was reported by Gierlinski et 
al. (1999), but they did not include the higher 
energy COMPTEL data. The inclusion of the 
COMPTEL data in the high state spectrum pro¬ 
vides evidence, for the first time, of a continuous 
power-law (with a photon spectral index of 2 . 6 ) 
extending beyond 1 MeV, up to ~ 10 MeV. No 
clear evidence for a cutoff in the power-law spec¬ 
trum can be discerned from these data. 

A power-law spectrum had also been observed 
by both OSSE and BATSE during the high X-ray 
state of February, 1994 (CGRO VP 318.1; Phlips 
et al. 1996; Ling et al. 1997). These earlier data 
correspond to the low level of hard X-ray flux near 
TJD 9400 in Figure 1. The spectrum observed 
during the 1994 high state showed a similar photon 
spectral index (F = 2.72 vs. T = 2.57 for the 1996 
high state spectrum), but the overall intensity of 
the power law was considerably lower (Gierlinski 
et al. 1999). Near 1 MeV, for example, the spec¬ 
tral amplitude was about 3 times lower in 1994 
than it was in 1996. This explains why Cygnus 
X-l was not observed by COMPTEL during the 
1994 high state. The extrapolation of the lower- 
intensity power-law fell below the sensitivity limit 
of COMPTEL. On the other hand, the intensity 
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observed in 1996 was sufficiently high to allow for 
a measurement of the spectrum by COMPTEL. 

We have used two different hybrid thermal/non- 
thermal Comptonization models (compps and 
eqpair) to fit broad band spectral data obtained 
in both the hard and soft spectral states. For the 
hard state analysis, we used data from CGRO cov¬ 
ering 20 keV up to 10 MeV. For the soft state anal¬ 
ysis, we augmented the CGRO data with lower 
energy data from BcppoSAX to provide improved 
constraints on the spectrum at energies down to 
0.5 keV. These fits provide a quantitative esti¬ 
mate of the electron distribution and allow us to 
probe the physical changes that take place dur¬ 
ing transitions between the low and high X-ray 
states. Hybrid Comptonization models have also 
been used to model the spectra of other black hole 
binaries in their soft state, such as GRS 1915+105 
(Zdziarski et al. 2001). 

The high energy spectrum of Cygnus X-l can¬ 
not be described by the bulk-motion Comptoniza¬ 
tion model alone, which predicts a sharp cutoff 
above ~ 100 keV (Laurent & Titarchuk 1999). 
The hybrid comptonization models provide an 
adequate fit to the data without requiring any 
contribution from bulk-motion Comptonization. 
Furthermore, the bulk-motion Comptonization 
power-law for L ~ 0.04Le corresponding to the 
soft state of Cygnus X-l (see below), was found 
by Laurent & Titarchuk (1999) to be very soft, 
with T ~ 3.5 at , i.e., much softer than the 
observed T ~ 2.5 (Gierlinski et al. 1999; Fron- 
tera et al. 2001). (See also the discussion in 
Zdziarski 2000.) Note that the XSPEC model 
of bulk-motion Comptonization, bmc (Shrader & 
Titarchuk 1998), does not include any high-energy 
cutoff and thus cannot be applied to our data (or 
any data extending to > 100 keV). 

Figure 9 shows a comparison of the spectra in 
the two states. For the hard state, we also show a 
typical spectrum at energies < 25 keV (BeppoSAX 
data from Di Salvo et al. 2001). We see that the 
two broad-band spectra cross each other at ~ 10 
keV and ~ 1 MeV. The dashed curve shows the 
model obtained by fitting the hard-state data from 
CGRO only (§3.1), and assuming Nu = 6 x 10 21 
cm -2 . We see that this model predicts the low- 
energy BeppoSAX data relatively well, underesti¬ 
mating somewhat the observed spectrum only at 
< 10 keV due to the presence of a pronounced 


soft X-ray excess present in the hard state (Ebi- 
sawa et al. 1996a; Frontera et al. 2001; Di Salvo et 
al. 2001), which is neglected in our model fitted to 
the data at > 20 keV. 

The bolometric flux or luminosity ratio between 
the soft state in June 1996 and the average for the 
hard state is ~ 4. This value is much more than 
the rough estimate of ~ 1.5 — 1.7 based on the 
ASM and BATSE occultation results (Zhang et al. 
1997), but is consistent with the results of Frontera 
et al. (2001), based on studies with BeppoSAX. 
Such a large value makes models of the state tran¬ 
sition based on a change of accretion rate plausi¬ 
ble. Given the larger luminosity in the soft state, 
the characteristic dimension of the hot plasma in 
the soft state based on the compactness fit is sim¬ 
ilar to that in the hard state, ~ 10 2 GM/c 2 . 

These data tend to support the general pic¬ 
ture that the transition between the hard and soft 
states results from a change in the disk transition 
radius between a hot inner corona (ADAF) and 
a cooler outer thin disk (e.g., Esin et al. 1998; 
Narayan et al. 1998; Poutanen 1998a,b; Poutanen 
& Coppi 1998). In the hard state, this transi¬ 
tion radius is relatively far from the black hole 
(at ~ 100 Schwarzschild radii). The spectrum 
is dominated by Comptonization off the thermal 
electrons in the hot inner corona. Radio emis¬ 
sion is also more pronounced in this state (Fender 
2001), with evidence for a radio emitting relativis¬ 
tic jet (Stirling et al. 2001). As the transition ra¬ 
dius moves inward, perhaps due to an increase in 
the accretion rate, the optically thick cool disk in¬ 
tercepts a larger fraction of the energy. The ther¬ 
mal energy dissipation in the corona is reduced 
considerably and the blackbody disk component 
(the principal component at soft X-ray energies) 
becomes more pronounced. 

Although our data tend to support the above 
picture, we have not attempted to model the ge¬ 
ometry in detail, since the precise geometrical con¬ 
figuration of the emitting region is largely un¬ 
known. Furthermore, our new data cover the en¬ 
ergy range near 1 MeV where geometry effects are 
difficult to study. One of the primary goals of the 
present paper is to determine the electron distri¬ 
bution of the radiating plasma. Our assumption 
of a spherical source geometry provides the neces¬ 
sary physics that is required to extract information 
on the electron spectrum. We have further pre- 
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sumed that the thermal and nonthermal electrons 
are in the same physical region. This assumption 
is based, in part, on the observations that show a 
negative correlation between the thermal and non¬ 
thermal components. This need not be the case 
in reality, however, but the present data cannot 
be used to determine the extent to which the two 
populations are co-located. A more detailed dis¬ 
cussion of geometrical effects in the context of the 
eqpair model, including Compton reflection and 
energy balance, can be found in Gierlinski et al. 
(1999). 

The shape of the electron distribution and its 
high energy tail can best be determined by mea¬ 
surements that extend into the MeV energy re¬ 
gion. The high-energy cutoff is related to the com¬ 
pactness of the source region, since it depends, in 
part, on the influence of 7 — 7 pair production. 
If 7 — 7 pair production is an important source 
of opacity, this would imply the presence of a sig¬ 
nificant level of pairs in the source region. In 
this way, a measure of the high-energy cutoff can 
help determine the nature of the emitting plasma 
(e — p or e ± ). Although a measure of e± annihi¬ 
lation radiation can also serve as a diagnostic of a 
pair plasma, it is likely that any annihilation radi¬ 
ation that may be present would be considerably 
broadened (and perhaps blue-shifted), and hence 
may not be readily observable. Measurements to 
date with HEAO-3 (Ling & Wheaton 1989) and 
with OSSE (Phlips et al. 1996) provide only up¬ 
per limits, or, at best, a marginal (1.9er) detection 
(Ling & Wheaton 1989) to the level of e ± anni¬ 
hilation radiation. This further underscores the 
need to define the high-energy cutoff as perhaps 
the best means for constraining the source com¬ 
pactness and the nature of the emitting region. 
If INTEGRAL, with its improved line sensitiv¬ 
ity, succeeds in measuring an annihilation feature, 
then constraints on the high-energy cutoff will be 
even more valuable. 

The presence of a non-thermal tail in the elec¬ 
tron distribution can also provide constraints on 
the strength of the magnetic field in the source re¬ 
gion. As pointed out by Wardzinski & Zdziarski 
( 2001 ), the presence of even a weak nonthermal 
electron tail increases strongly the emissivity of 
the cyclo-synchrotron process with respect to the 
pure thermal case. If the Compton-scattering elec¬ 
trons in Cygnus X-l were purely thermal, that 


process appears in general to be too inefficient to 
provide all of seed photons for the Comptoniza- 
tion under simple assumptions of equipartition 
(Wardzinski & Zdziarski 2000). Since we do see a 
blackbody component at low energies (Ebisawa et 
al. 1996a; Di Salvo et al. 2001), this inefficiency is 
consistent with the seed photons for Comptoniza- 
tion provided by the blackbody rather than by the 
cyclo-synchrotron photons. On the other hand, 
the tail parameters obtained by McConnell et al. 
( 2000 a) yielded such a copious supply of cyclo¬ 
synchrotron seed photons that the corresponding 
luminosity would become ~ 10 2 times that ob¬ 
served (Wardzinski & Zdziarski 2001). This con¬ 
clusion is confirmed for the tail parameters fit¬ 
ted here (G. Wardzinski, private communication). 
Thus, either the magnetic field in Cygnus X-l is 
substantially below equipartition (at least an or¬ 
der of magnitude) or the observed photon tail has 
a different origin than that due to a high energy 
electrons. In either case, this has important impli¬ 
cations for models of the accretion flow in Cygnus 
X-l. 

These studies also have implications that go 
beyond that of studying individual black hole 
sources. Given the close spectral similarity be¬ 
tween black-hole binaries in the hard state and 
Seyferts (e.g., Zdziarski 2000), it is possible that 
similar tails are present in the spectra of the latter 
objects. Stecker et al. (1999) have suggested that 
the hard tail emission seen in sources like Cygnus 
X-l might account for an important component 
of the cosmic diffuse background radiation in the 
200 keV - 3 MeV energy band (see also Stecker 
2001). Note, however, that the tail of Cygnus 
X-l above 1 MeV contains relatively little flux, 
1.3% of the bolometric (model) flux, for the fit 
with eqpair. If a similar value is characteristic 
of Seyferts, the combined emission from their high 
energy tails may be too weak to account for the ob¬ 
served extragalactic MeV background, perhaps ar¬ 
guing against the proposal by Stecker et al. (1999). 

The next major satellite for this energy range, 
INTEGRAL, is expected to have only slightly bet¬ 
ter continuum sensitivity than COMPTEL at en¬ 
ergies near 1 MeV with both its IBIS and SPI ex¬ 
periments (Schonfelcler et al. 2001). Furthermore, 
the much narrower FoV of the INTEGRAL instru¬ 
ments (~ 15°) will mean that there will likely be 
only a limited number of observations of Cygnus 
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X-l. This is in stark contrast to the COMPTEL 
situation, in which the large FoV of COMPTEL 
(~ 60°) resulted in many weeks of exposure, most 
of which were obtained during the low X-ray state. 
Given the large low-state exposure of CGRO , it is 
quite likely that INTEGRAL may not be able to 
offer any significant improvement in our knowledge 
of the hard state continuum spectrum at MeV en¬ 
ergies. The CGRO data may therefore provide 
the best view of the hard state MeV continuum 
for many years to come. However, COMPTEL 
is very limited in the data that it collected for 
the soft state spectrum. Additional soft state ob¬ 
servations with INTEGRAL could therefore prove 
valuable. An important goal would be to search 
for a cutoff in the energy spectrum. Pinning down 
the energy of this cutoff would be a very impor¬ 
tant next step in our understanding of the high 
energy spectrum of Cygnus X-l. In this regard, 
INTEGRAL may be an extremely useful tool for 
collecting additional soft state spectral data, pro¬ 
viding that suitable target-of-opportunity obser¬ 
vations can be acquired. 
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Fig. 1. - X-ray time histories of Cygnus X-l cov¬ 
ering nearly the entire CGRO mission. The hard 
X-ray data come from BATSE data that are de¬ 
rived from Earth occultation analysis in the 20- 
100 keV energy range. The soft X-ray data (2-10 
keV) are from the All-Sky Monitor (ASM) on the 
Rossi X-Ray Timing Explorer (RXTE). 



Truncated Julian Day (TJD) 

Fig. 2.— The time interval of CGRO Viewing 
Period 522.5 is shown relative to the hard (upper) 
and soft (lower) X-ray time histories from BATSE 
and the RXTE-ASM, respectively. Note the very 
rapid transition into and out of the high X-ray 
state as seen in soft X-rays. 
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Fig. 3.— COMPTEL imaging of the Cygnus re¬ 
gion as derived from 1 -3 MeV data collected dur¬ 
ing the soft state of June, 1996 (CGRO viewing 
period 522.5). The outer contour (dotted line) 
represents the effective FoV of COMPTEL (with 
a 30° radius). The remaining contours represent 
constant values of the quantity —2 In A, where A 
is the likelihood ratio. The contours start at a 
value of 15, with a step size of 5. The likelihood 
reaches a value of 30.1 at the location of Cygnus 
X-l (denoted by the diamond). 
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Fig. 4.— The 1, 2 and 3 -a location contours de¬ 
rived from the likelihood map in Figure 3. The 
emission is consistent with a point source at the 
location of Cyg X-l, with no significant contribu¬ 
tion from PSR 1951+32, a pulsar that has been 
detected in a timing anlaysis of COMPTEL data 
(Kuiper et al. 1998). 



E [keV] 

Fig. 5.— The average CGRO spectrum of Cygnus 
X-l in the hard state fitted with the eqpair model 
(solid curve). Data points from BATSE and OSSE 
are represented as blue open circles and red aster¬ 
isks, respectively. COMPTEL data are shown as 
thick crosses. All the data are normalized to that 
of OSSE. Upper limits have been removed for the 
sake of clarity. 
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Fig. 6.— Components of the eqpair fit for the 
hard state. All spectra are intrinsic, i.e., corrected 
for absorption. The long dashes, short dashes, 
and dots correspond to the unscattered blackbody, 
scattering by thermal electrons, and Compton re¬ 
flection, respectively. The solid curve is the to¬ 
tal spectrum. Scattering by the nonthermal elec¬ 
trons accounts for the high-energy tail above the 
thermal-Compton spectrum given by the short 
dashes, starting at ~ 1 MeV. 
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Fig. 7.— The simultaneous BeppoSAX-CGRO 
spectrum of Cygnus X-l in the soft state fitted 
with the eqpair model (solid curve). Included are 
data from the LECS (green), HPGSPC (cyan open 
squares) and PDS (yellow open triangles) instru¬ 
ments on board BeppoSAX and from the OSSE 
(red asterisks), BATSE (blue open circles) and 
COMPTEL instruments on CGRO. All the data 
are normalized to that of OSSE. 
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Fig. 9.— Comparison of the spectra in the hard 
and soft state of Cygnus X-l, as fitted with the 
eqpair model (solid curves). All data are nor¬ 
malized to that of OSSE. The hard-state data be¬ 
low 25 keV represent a typical hard-state spectrum 
(the BeppoSAX data of Di Salvo et al. 2001), and 
the dashed curve shows the corresponding model 
obtained by fitting the CGRO data only. 


Fig. 8.— Components of the eqpair fit for the 
soft state. All spectra are intrinsic, i.e., corrected 
for absorption. The long dashes, short dashes, 
dot/dashes and dots correspond to the unscattered 
blackbody, scattering by thermal electrons, the 
scattering by nonthermal electrons, and Compton 
reflection/Fe Ka fluorescence, respectively. The 
solid curve is the total spectrum. 
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Table 1 

COMPTEL Flux Measurements for VP 522.5. 


Energy 

(MeV) 

Counts in 
Dataspace 

Source 

Counts 

Flux 

(cm -2 s -1 MeV -1 ) 

0.75 - 1.0 

33,066 

177 ± 202 

5.4(±6.1) x 10 -4 

1.0 - 1.4 

88,086 

697 ± 258 

5.3(±2.0) x 10 -4 

1.4 - 2.0 

107,010 

756 ± 265 

2.6(±1.0) x 10 -4 

2.0 - 5.0 

157,457 

1092 ±273 

5.8(±1.5) x 10 -5 

5.0 10.0 

30,158 

191 ±90 

5.8(±2.7) x 10 -6 

10.0 - 30.0 

7,967 

31 ±31 

3.2(±3.2) x 10 -7 
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Table 2 

Parameters of the hybrid models for the hard and soft state 


N H 

kT s 

4 

en/e. 

^nth/^h 

7 min 

P , r inj kT 

Ti 

T 

0/27T 

t 

F 

L 

x 2 /- 


hard state, compps 

6f 

0.2f 

- 

- 

- 

1 QQ + 051 

l.oy_ 0 34 

5.4ig;| 58^22 

- 

9 q+ 0-6 

z<y -0.4 

0.5f 

- 

3.38 

1.62 

37/48 

hard state, eqpair 

6f 

0.2f 

1 oT2.5 

1 + 3 

0.082!°;°“ 

1.5f 

2.0^2 0 9 ° a 

1 34+ 0 - 40 
J--Oft_ 0 .5o 

1.45 a 

0 52+ 0 - 06 

U.OZ_Q Q 5 

- 

3.56 

1.70 

31/46 

soft state, compps 

6 0" 1 " 0 ' 2 
D,u -0.1 

0 39" 1 " 0 ' 01 
u,oy _0.01 

- 

- 

- 

1 qq+0.12 

1.03_ 0 12 

3 -5±S:l 6311 

- 

n 10+0-04 
U-io_Q 02 

1 4+0-3 

-‘-•^-0.5 

290" 1 " 350 
^ —180 

13.1 

6.3 

199/239 

soft state, eqpair 

o o 
+ 1 
o 
cd 

0 q7+ 0 - 01 
u " 3 ' -0.01 

o o+38 
J.Z-2 i 

o ]7+ 0 - 01 
U -- L * -0.01 

0.68±°; 2 i° 

1.5f 

2.6111 65 a 

0 11 + 0.02 
uxx - 0.02 

0.11 a 

1 q+0.3 

1 -0.3 

100 +21 ° 
xuu —60 

13.2 

6.3 

199/238 


Note.—Ah, kT, and £ are in units of 10 21 cm -2 , keV, and ergcms -1 , respectively. F , the unabsorbed model bolometric flux using the normalization of the OSSE 
spectrum, and L, the corresponding luminosity assuming isotropy and a distance of 2 kpc are in units of 10 —8 erg cm -2 s —1 and 1() 37 ergs -1 , respectively. Parameters fixed 
in the fit are denoted by ‘f\ The single-parameter uncertainties correspond to a 90% confidence, i.e., A% 2 = 2.71. 

a The electron temperature and total optical depth calculated from energy and pair balance for the best-fit model (i.e., not free parameters). 




